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Abstract. We study the spreading of correlations and other physical quantities 
in quantum lattice models with interactions or hopping decaying like r~°‘ with the 
distance r. Our focus is on exponents a between 0 and 6, where the interplay of 
long- and short-range features gives rise to a complex phenomenology and interesting 
physical effects, and which is also the relevant range for experimental realizations 
with cold atoms, ions, or molecules. We present analytical and numerical results, 
providing a comprehensive picture of spatio-temporal propagation. Lieb-Robinson- 
type bounds are extended to strongly long-range interactions where a is smaller than 
the lattice dimension, and we report particularly sharp bounds that are capable of 
reproducing regimes with soundcone as well as supersonic dynamics. Complementary 
lower bounds prove that faster-than-soundcone propagation occurs for a < 2 in any 
spatial dimension, although cone-like features are shown to also occur in that regime. 
Our results provide guidance for optimizing experimental efforts to harness long-range 
interactions in a variety of quantum information and signaling tasks. 


1. Introduction 

Traditionally, the study of lattice models has focused on Hamiltonians where interactions 
and/or hopping is restricted to a few neighboring sites. Only recently there has been a 
surge of interest in long-range interacting systems where interaction strengths or hopping 
amplitudes decay like a power law at large distances r. This interest was triggered 
on the experimental side by progress in the control of ultra-cold atoms, molecules, and 
ions, which led to the realization of a variety of long-range systems. Examples include 
magnetic atoms [1], polar molecules [2], trapped ions [3, 4, 5, 6], Rydberg atoms [7], 
and others. On the theoretical side, intriguing physical effects and properties have 
been predicted for long-range interacting quantum systems, including nonequivalent 
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statistical ensembles and negative response functions [8, 9], equilibration time scales 
that diverge with system size [10, 11, 12], prethermalization [13, 14], and others. 

In this article we study the propagation in time and space of various physical 
quantities, and this is another topic where long-range interactions lead to peculiar 
behavior. A number of papers devoted to this topic have appeared in the past two 
years, reporting results on the spreading of correlations, information, or entanglement 
[15, 16, 17, 18, 19, 20, 21]. In short-range systems, all these quantities are known to 
propagate approximately within a soundcone, reminiscent of the lightcone in relativistic 
theories, with only exponentially small effects outside the cone. This behavior is termed 
quasilocality and was rigorously proved by Lieb and Robinson for a class of short-range 
interacting lattice models [22] . In the presence of long-range interactions this picture is 
altered signihcantly: the concept of a group velocity breaks down, and the spreading of 
correlations, information, or entanglement may speed up dramatically. This, in turn, 
has a bearing on all kinds of dynamical properties, and one might hope to harness long- 
range interactions for fast information transmission, improved quantum state transfer, 
or other applications. 

Much of our understanding of propagation in long-range systems comes from 
analytical or numerical studies of model systems, where for example correlations or 
entanglement between lattice sites i and j are calculated as functions of time t and 
spatial separation d{i,j). Typical examples of such results, similar to some of those in 
[15, 16, 17, 18, 19, 20, 21], are shown in hgure 1 for a number of different models, physical 
quantities, and exponents a. For larger a (hgure 1 right), the behavior is reminiscent of 
the short-range case, with only small effects outside a cone-shaped region. For small a 
(hgure 1 left), correlations propagate faster than any hnite group velocity would permit, 
and are mostly conhned to a region with power law-shaped boundaries. For intermediate 
a (hgure 1 center), a crossover from cone-like to faster-than-cone behavior is observed. 
While these three regimes seem to be typical and occur in many of the models studied, 
notable exceptions (some of which will be discussed further below) do occur and lead 
to a more complicated overall picture. 

Besides model calculations, Lieb-Robinson-type bounds have contributed signih¬ 
cantly to our understanding of propagation in long-range interacting models. The hrst 
result of this kind. 

Ml \B\ _ i) 

\\[OA{t),OBm\\ < CIIO^II IIObII (1) 

valid for exponents a larger than the lattice dimension D, was reported by Hastings 
and Koma [23]. Here, A,BgA are non-overlapping regions of the lattice A, and 
Oa{0) and Ob( 0) are observables supported only on the subspaces of the Hilbert space 
corresponding to A and B, respectively. H-H denotes the operator norm, and d{A,B) is 
the graph-theoretic distance between A and B ]:. The relevance of the bound (1) lies in 

f The graph-theoretic distance is the number of edges along the shortest path connecting the two 
regions. 
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Figure 1. Propagation patterns as a function of distance 6 = d{i,j) and time t for 
different long-range exponents a. To highlight the generality of the phenomena we 
discuss in this article, we use different models and physical quantities as examples. 
Left: For a long-range Ising chain with a = 1.2, we show the probability to detect 
a signal sent through a quantum channel from site 0 to 5 [15]. The green line is a 
guide to the eye and shows a power law d oc P''’. Center: Connected equal-time 
correlations between lattice sites 0 and 5 in a long-range field theory in one spatial 
dimension with a = 4 [21]. After an initial cone-like spreading, a cross-over to power 
law-shaped contours is observed. The green dashed curve is a guide to the eye. Right: 
The spreading of entanglement as captured by the mutual information between two 
lattice sites separated by a distance 5 in the long-range hopping model (13) with a = 8, 
starting from a staggered initial state (see text). Entanglement is sharply confined to 
the interior of a cone. 


the fact that a number of physically interesting quantities, like equal-time correlation 
functions, can be related to the operator norm of the commutator on the left-hand 
side of (1), so that similar bounds hold also for these physical quantities [24, 25]. For 
any a, a contour plot of the bound (1) looks qualitatively like the plot in Fig. 1 (left), 
although with logarithmic contour lines instead of power laws. This implies that, while 
correct as a bound for all a > D, the shape of the propagation front (hgure 1 center 
and right) is not correctly reproduced by (1) for intermediate or large values of a. 
Another bound put forward in [26] improves the situation for the case of large a, but 
turns out to be weaker than (1) for smaller values §. Summarizing the situation, the 
existing Lieb-Robinson-type bounds struggle to reproduce the transition from cone-like 
to faster-than-cone propagation for intermediate a as in hgure 1 (center) ||. For small 
a, no bounds have been published so far. 

In this article we prove general bounds, complemented by model calculations, 
that help to establish a comprehensive and consistent picture of the various kinds of 
propagation behavior that occur in long-range interacting lattice models. We extend 
Lieb-Robinson-type bounds to strong long-range interactions where a < D. This is 
complemented by model calculations showing that, even in the regime a < D of strong 
long-range interactions, cone-like propagation may be a dominant feature. We also prove 
that faster-than-cone propagation can occur for all a < 2 in any spatial dimension, and 

§ See Appendix A.3 for a more detailed discussion of the bound in [26]. 

II We could not compare the tightness of the matrix exponential bound with that of the bound in [27], 
as several of the constants occurring in that bound were not specified. 
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this answers a question put forward in [6]. For intermediate exponents a, we advocate 
the use of a Lieb-Robinson-type bound in the form of a matrix exponential, which is 
tight enough to capture the transition from a cone-like to a faster-than-cone propagation 
as in hgure 1 (center), and is also computationally efficient. 


2. Lieb-Robinson bounds for a. < D 


For deriving analytical results in the regime a < D, an understanding of the time scales 
of the dynamics turns out to be crucial. The presence of strong long-range interactions 
is known in many cases to cause a scaling of the relevant time scales with system 
size [28, 10, 11, 12, 13, 14]. For long-range quantum lattice models the fastest time 
scale oc N~'^ was found to shrink like a power law with increasing system size N, 
where g is a positive exponent [12, 13]. This observation makes clear why previous 
attempts to derive a Lieb-Robinson-type bound for a < D failed: in the large-limit 
the dynamics becomes increasingly faster, and hence propagation is not bounded by any 
hnite quantity. Considering evolution in rescaled time r = tN'^ can resolve this problem 
and allows us to obtain a hnite bound in the thermodynamic limit. 

On an arbitrary H-dimensional lattice A with N sites we consider the Hilbert space 

N 

Ar = (g)^ (2) 

i=l 

with hnite-dimensional local Hilbert spaces On a generic Hamiltonian 

W = E (3) 

XcA 

with n-body interactions is dehned, with local Hamiltonian terms hx compactly 
supported on the hnite subsets X C A. The Hamiltonian is required to satisfy the 
following two conditions. 


Boundedness, 


- 

XBi,j 


A 

[1 + d{i,j)]^ 


( 4 ) 


with a hnite constant A > 0. This condition, also used in [23], is a generalization 
of the dehnition of power law-decaying interactions, and it reduces to the usual 
dehnition in the case of pair interactions, i.e., when X consists only of the two 
elements i and j. 


(ii) Reproducibility, 

‘VaE 


< 


P 


fceA 


[1 + d{i,k)Y[l + d{k,j)Y [1 + d{i,j)Y 


( 5 ) 


for hnite p > 0, with 

aAa = 1/ sup 




( 6 ) 
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Figure 2. Bound (8) in physical (not rescaled) time t for a = 1/2 and lattice sizes 
N = 10^, 10^, and 10^ (from left to right), illustrating the speed-up of the propagation 
with increasing lattice size. For simplicity all constants in (8) are set to unity. 


The lattice-dependent factor is the same that is frequently used to make a long-range 
Hamiltonian extensive [29, 10], but we use it here for a different purpose. Asymptotically 
for large regular lattices, one hnds [ 10 ] 

( for 0 ^ a < A), 

-Aa ~ < C 2 /\riN for a = H, (7) 

C 3 for a > D, 

with a-dependent positive constants ci, C 2 , and C 3 . Eq. (5) is a modihed version of one 
of the requirements for the proof in [23], but due to the modihcation by the factor .Aa 
the condition is satished for a larger class of models, including regular H-dimensional 
lattices with power law-decaying interactions with arbitrary positive exponents a [30]. 
For the above described setting we derive in Appendix A.2 the Lieb-Robinson-type 
bound 

\ A \ I Rl (pAA — 1 i 

11 [OA(r^A), Ob(0)] 11 < C I |Oa| I I |Qb 11 ^ (8) 

in rescaled time 

T = tjPpK- (9) 

This bound reproduces qualitative features of supersonic propagation (as in hgure 1 
left), and also accounts for the system-size dependence of the time scale of propagation 
for exponents a < D. While the bound ensures well-dehned dynamics in rescaled time r 
in the thermodynamic limit, it describes a speed-up in physical time t of the propagation 
with increasing lattice size, as illustrated in hgure 2 . 

3. Matrix exponential bounds for intermediate ck 

For long-range models with intermediate exponents, in the range 3 < a < 6 or even a 
bit larger, one observes an interplay of cone-like and supersonic propagation (hgure 1 
center). This is the most relevant regime for experimental realizations of long-range 
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Figure 3. Spacetime plots of the matrix exponential bound (11) for several values 
of a in a one dimensional system with L = 201 lattice sites and periodic boundary 
conditions. Left: for a = 1.2 the bound recovers a propagation front with a shape 
similar to the one of the Ising model in figure 1 (left). Center: for intermediate a = 4 
a transition from soundcone to supersonic dynamics is being heralded. Right: the two 
regimes of soundcone-like and supersonic dynamics are fully exposed for a = 8. 


interactions by means of cold atoms or molecules, but a theoretical description of 
the shape of the propagation front turns out to be challenging. Existing bounds [26] 
are discussed in Appendix A.3. Here we report bounds that capture the features of 
the propagation front as observed in long-range models with intermediate exponents, 
showing a clear and sharp crossover from cone-like to supersonic propagation. 

As in section 2, our setting is a D-dimensional lattice A consisting of N sites and 
a Hilbert space (2) with hnite-dimensional local Hilbert spaces. We consider a generic 
Hamiltonian with pair interactions, 

1 ^ 

^ k^i 

where the pair interactions hki are bounded operators supported on lattice sites k and I 
only. As observables Oa and Ob consider bounded operators that are supported on 
single sites A = {i} and B = {j}. In this setting, we prove in Appendix A.l a bound 
in the form of an Af x Af matrix exponential, 

||[Oi(t),Oj(0)]|| < 2||Oi|| IIOjII (^exp[2KJ \t\]. - - , (11) 

where J is the interaction matrix with elements 

Jk,i = \\hki\\ (12) 

and n = sup^gy^ Yhk ^ri,k- In one-dimensional homogeneous lattice models the interaction 
matrix J is of Toeplitz type and thus (11) can be evaluated in time using the 

Levinson algorithm [31]. For translationally invariant one-dimensional systems, J is 
a circulant matrix, which permits an analytical solution of (11) by means of Fourier 
transformation. 

The bound (11) is tighter than the bounds in [23, 25, 26], and the crossover from 
cone-like to supersonic propagation is nicely captured (see figure 3). Due to its form as a 
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matrix exponential, the bound is less explicit than others, and asymptotic properties are 
not easily read off. But since the calculation of a matrix exponential scales polynomially 
in the matrix dimension N (like or even faster [32]) the bound can easily be 

evaluated for large lattices up to on a desktop computer. This is orders of 

magnitude larger than the sizes that can be treated by exact diagonalization, and covers 
the system sizes that can be reached for example with state-of-the-art ion trap based 
quantum simulators of spin systems [3]. Different from other bounds of Lieb-Robinson- 
type, our matrix exponential bound is computed for the exact type of interaction matrix 
realized in a specihc experimental setup. This improves the sharpness of the bound, and 
can make it a useful tool for investigating all kinds of propagation phenomena in lattice 
models of intermediate system size. 

4. Long-range hopping for small ck 

The bounds discussed in sections 2 and 3 are valid for arbitrary initial states, and 
therefore it may well happen that propagation for a given model and some, or even 
most, initial states is signihcantly slower than what the bound suggests. Indeed, linear 
(cone-like) propagation was observed in model calculations even for moderately large 
exponents like a = 3 [17, 16, 18, 15, 19]. But, as we show in the following, such cone-like 
propagation can, for suitably chosen initial states, even persist into the strongly long- 
range regime 0 < a < D. In this and the next section we analyze free fermions on a 
one-dimensional lattice with long-range hopping, which is arguably the simplest model 
to illustrate cone-like propagation in long-range models and to explain the observation 
on the basis of dispersion relations and density of states. While strictly speaking such 
a long-range hopping model does not meet the conditions under which Lieb-Robinson 
bounds have been proved, it proves helpful for understanding the conditions under 
which cone-like propagation may or may not be observed in other long-range interacting 
models. 


4.1. Long-range hopping model 


Consider a free fermionic hopping model in one dimension with periodic boundary 
conditions, 

N N-1 

^ = “2 ’ (^ 3 ) 

i=l 1=1 


where cj, Cj are fermionic creation and annihilation operators at site j. We choose 
long-range hopping rates oc df^, where 


di = 


I 


N-l 


if / < N/2, 
if I > N/2, 


(14) 
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is the shortest distance between two sites on a chain with periodic bonndary conditions. 
A Fonrier transformation brings the Hamiltonian into diagonal form 



(15) 


k 


with 



(16) 


and dispersion relation 



( 17 ) 


where k = 2TTm/N with m = 1 ,..., A^. 

4-2. Propagation from staggered initial state 

We choose a staggered initial state |1010 ...) in position space, i.e., initially every second 
site is occupied. For simplicity of notation we assume the number N of lattice sites to 
be even. A straightforward calculation, similar to that in [33], yields 



(18) 


n=l 

for the time-dependence of the occupation number at lattice site j, where 



(19) 


and k = 2'Km/N with m = l,...,iV. In figure 4 (left) the time evolution of the 


occupation number {nj{t)) is plotted for different values of a, showing that the time it 
takes to relax to the equilibrium value of 1/2 increases dramatically for small a (note the 


logarithmic timescale). This may seem counterintuitive, as a longer interaction range 


may naively be expected to lead to faster propagation. The effect can be understood 


from figure 4 (right), showing the spectrum of the frequencies A in the cosine terms of 
Eq. (18). As a decreases, the majority of these frequencies lie within a small window 
around zero, implying very slow dephasing of the cosine terms. 

A more refined picture of the propagation behavior can be obtained by studying the 
spreading of correlations. Starting again from a staggered initial state, a straightforward 
calculation similar to that in [34, 33], and similar to the one leading to (18), yields 



( 20 ) 


k 


Figure 5 (bottom) shows contour plots in the (5, t)-plane of the absolute values of the 
correlations (20) for different values of a. For all a shown, a cone-like propagation front 
is clearly visible, even in the case of a = 3/4 < H. Two properties of the cone can be 
observed to change upon variation of the exponent a\ (i) The boundary of the cone is 
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{^j) 



A 



Figure 4. Left: Time dependence of the occupation number of site j for different a, 
starting from a staggered initial state. Right: A as a function of k. The system size is 
N = 200 in both plots. 



Figure 5. Contour plots in the (<5, t)-plane, showing correlations (20) between sites 0 
and 5 in the fermionic long-range hopping model for N = 200 lattice sites and various 
values of a, starting from a staggered initial state. 


rather sharp for larger a (like a = 3), whereas correlations “leak” into the exterior of 
the cone for smaller a (like a = 3/4). (ii) The velocity of propagation, corresponding to 
the inverse slope of the cone, decreases with decreasing a [see hgure 6 (left)], conhrming 
the counterintuitive observations of hgure 4 (left). We will argue in section 4.4 that 
some of these features can be understood on the basis of the dispersion relation (17) 
and the density of states of the long-range hopping model. 

4-3. Dispersion and group velocity 

In the limit of large system size the dispersion relation takes the form 

e(fc) = -[Li,(e''=)+Li„(e-''^)] (21) 

where Lie is the polylogarithm [35], and this function is plotted in hgure 7 (left). For 
a = 3 the dispersion e is a smooth function of k, while it shows a cusp at A: = 0 for 
a = 2, and a divergence at fc = 0 for a = 1. Correspondingly, the derivative e'{k) as 
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shown in fignre 7 (right) is discontinnous at fc = 0 for a = 2, and diverges at /c = 0 
for a = 1. More generally we can analyze e' in the vicinity of fc = 0 by considering the 
difference qnotient between the zeroth and the first mode, 


e{27i/N) - e(0) 

N ^ \cos{27il/N) - 1| 

27r(l - 0)/iV 

27r df 

1=1 ' 


“ Sir df iV ^ 


( 22 ) 


(23) 


In the large-iV limit we approximate the snm by an integral, 

A pN/2 c\ 

nL = ( 24 ) 

This implies that, for a < 2, the derivative e' diverges at /c = 0 in the limit of infinite 
system size. Interpreting e'(0) as a gronp velocity, we infer that we have a finite gronp 
velocity only for a > 2, whereas the concept of a group velocity breaks down for a < 2 
This finding can help us to understand figure 5: For a > 2 a finite group velocity restricts 
the propagation to the interior of a cone, which makes this cone appear rather sharp. For 
a < 2, although a cone is still visible, larger (and, in fact, arbitrarily large) propagation 
velocities may occur and are responsible for the “leaking” of correlations outside the 
cone. 

The threshold value a = 2 for supersonic propagation (i.e., propagation not 
bounded by any finite group velocity) is also found in a different context, by very different 
methods. In [15] it was proved that information can be transferred supersonically 
through a quantum channel with hnite local Hilbert space dimension for any a < D + 1, 


^ The same conclusions about dispersion relations and group velocities also hold for long-range 
interacting XX and XXZ spin models when restricting the dynamics to the single magnon sector, 
as the dispersion relations of these models are essentially identical to (17). 


P 




Figure 6. Left: Dominant velocity of propagation, as read off from the inverse slope of 
the striking cones in figure 5, plotted as a function of the exponent a. Right: Density 
of states (25) for a = 1, 2 and oo. 
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Figure 7. Dispersion relation (21) (left) and its derivative e'(fc) (right) for the long- 
range fermionic hopping model (13) with exponents a = 1, 2, and 3. 


while no such proof exists for a > D + 1 but this result requires the measuring 
of observables supported on semi-infinite sublattices, which is not the most physical 
scenario. In Appendix B we use techniques similar to those in [15], but apply them to a 
model with translationally invariant interactions, to prove that supersonic transmission 
through a quantum channel can occur for any a < 2, also for measurements performed 
on single lattice sites. This result can be seen as complementary to the experimental 
observations in [6], where supersonic propagation of correlations was observed for 
exponents up to a ~ 1.19 in a one-dimensional lattice. 

4 . 4 . Density of states and typical propagation velocities 

From figure 5 and the discussion in section 4.3 we have seen that, while supersonic 
propagation can occur for a < 2, cone-like propagation is observed for these values of a 
at least for some initial states. In this section we will argue that the qualitative features 
of the observed behavior can be understood on the basis of the density of states 

= ^is)" 

in the large system limit. Eq. (25) can be rewritten as 

dk (26) 

where the sum is taken over all roots of the argument of the delta function. The 
polylogarithms that appear in the dispersion relation (21) can be analytically evaluated 




Hk-ko) 




■*" For models with infinite dimensional local Hilbert spaces JPi, supersonic propagation can occur also 
in models with nearest-neighbor interactions, although this appears to happen only under rather specific 
circumstances [36]. 
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for certain integer valnes of a, yielding 


/ 1 


for a = 1, 



for a = 2, 


(27) 


1 


for a —)■ oo, 


^ \/4: — 


where 0 is the Heaviside step fnnction. For those three valnes of a, the density of states 


is plotted in hgnre 6 (right), bnt other cases can be evalnated nnmerically (not shown 
in the hgnre). Again, as for the gronp velocity in hgnre 7 and the classical information 
capacity in Appendix B, we hnd a threshold valne of a = 2, as explained in the following. 

For a < 2, the density of states p is nonzero for all n, implying that propagation is 
not bonnded by any hnite maximnm velocity. The maximnm of p, however, is at n = 0 
for all a < 2, and this gives an indication that slow propagation with a small velocity 
is favored, althongh larger velocities do occnr [as in hgnre 5 (left and center)]. The 
maximnm at n = 0 becomes more sharply peaked when a approached zero, explaining 
the vanishing of the inverse slope of the cone in hgnre 5 in that limit, as shown in hgnre 6 


(left). 


For a > 2, the density of states p is nonzero only on a hnite interval [—Wmax, T'^max], 
where n^ax depends on a. For a > 2 the density of states diverges, and therefore takes 
on its maximnm, at iwmax- This implies that the maximnm velocity is favored, althongh 
smaller velocities also occnr [as in hgnre 5 (right)]. 

5. Conclusions 

In this paper we have stndied, from several diherent perspectives, the noneqnilibrinm 
dynamics of lattice models with long-range interactions or long-range hopping, and 
in particnlar the propagation in space and time of correlations and other physical 
qnantities. The focns of onr work is on the competition between linear, cone- 
like propagation and faster-than-linear, snpersonic propagation. We illnstrate this 
competition in two regimes, both relevant for experimental realizations of long-range 
many-body systems in cold atoms, ions, or molecnles: 

(i) For small exponents a < 2 we prove that snpersonic propagation can occnr. At 
the same time, in snch systems cone-like spreading can be the dominant form 
of propagation, with snpersonic ehects appearing only as small corrections [as in 
hgnre 5 (center)]. 

(ii) For intermediate exponents (ronghly between 3 and 8), propagation is observed to 
be linear initially, with snpersonic ehects setting in at larger times and distances 
[as in hgnre 1 (center) and hgnre 3 (right)]. 

To explain these observations, we provide model calcnlations as well as general bonnds 
that provide a comprehensive and consistent pictnre of the varions shapes of propagation 
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fronts that can occnr. Two of onr resnlts are Lieb-Robinson-type bonnds, valid for 
large classes of models with long-range interactions. The first is a bonnd for models 
with exponents a smaller than the lattice dimension D, a regime for which hitherto 
no such bounds existed. Key to deriving the bound is the insight that for a < D the 
propagation speed in general scales asymptotically like a power law with the system 
size, and a meaningful bound therefore has to be derived in rescaled time r as defined in 
(9). In physical time t, the bound then describes the increase of the propagation speed 
with increasing lattice size, as illustrated in figure 2. The second Lieb-Robinson-type 
bound we report is essentially a cheat, as we stop half way through the derivation of a 
“conventional” Lieb-Robinson bound. Specializing this result to single-site observables 
and Hamiltonians with pair interactions only, we obtain an expression that can be 
evaluated numerically in an efficient way, easily reaching system sizes of 1^(10'^). This 
bound (11) is sharp enough to capture cone-like as well as supersonic behavior. In 
experimental studies of propagation in long-range interacting lattice models [6, 5], the 
currently feasible lattice sizes are small and measured data can be compared to results 
from exact diagonalization. However, experimental work on systems of larger size is in 
progress, and exact diagonalization will not be feasible in that case. We expect that the 
matrix exponential bound (11) can provide guidance and sanity checks when analyzing 
the results of such experiments. 

In the second half of the paper we complemented the bounds with results of one 
of the simplest long-range quantum models, namely a fermionic long-range hopping 
model in one dimension. We observed that cone-like propagation fronts can be a 
dominant feature also for small values of a, and we explain the opening angle of such 
a cone, as well as the interplay of cone-like and supersonic features, on the basis of the 
dispersion relation combined with the density of states. These results indicate that it will 
depend crucially on the fc-modes occupied whether cone-like or supersonic propagation 
is dominant. We expect that such an improved understanding can provide guidance 
for optimizing experimental efforts to harness long-range interactions in a variety of 
quantum information and signaling tasks. 


Appendix A. Lieb-Robinson bounds 


Appendix A.l. Derivation of the matrix exponential bound 


As in section 2, we consider a D-dimensional lattice A consisting of N sites, a Hilbert 
space (2) consisting of hnite-dimensional local Hilbert spaces, and a generic Hamiltonian 
with pair interactions (10). Let Oa and be two bounded linear operators compactly 
supported on subsets A,B(ZA with A fl D = 0. Under these conditions, similar to the 
derivation of Eq. (2.12) of Ref. [25], one can derive the upper bound 


|||O.4W.Ob(0)]||<2||O^||||Ob||5^ 

71 = 1 


(2|t|) 


(A.l) 


nl 
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For pair interactions, and considering observables Oa and Ob supported on single lattice 
sites only (i.e., A = {z} and B = {j}), the coefficients a„ are upper bounded by 

On ^ ^ Ji,ki'^ki,k2 ■ ■ ■ Jk„-i,j = K {J , (A-2) 

/Cl . jkfi—i 


where J is the interaction matrix with elements Jk,i = \\hki\\ and k = 
Then (A.l) can be written as 


||[Oa(^), Qb(0)]|| ^ (2^ 1^1) 

2 ||Oa|| \\Ob\\ “ V^ n\ 


exp {2Jn\t\)..-Sij, 


(A.3) 


which proves (11). 

For translationally invariant one-dimensional lattices, J is a circulant matrix and 
can be diagonalized by means of Fourier transformation. For the example of power law 
interactions hki oc d{k,l)~°', the diagonal elements of the Fourier-transformed matrix J 
are given by 


e{k) 


L{iv- 1 )/ 2 J 


cosink) Q^kN /2 


n=l 


n“ 2{N/2Y 

if N even 


(A.4) 


with k = 27rm/N, 0 < m < N. Using these eigenvalues, J can be exponentiated in 
the diagonal basis, followed by an inverse Fourier transformation to evaluate the Lieb- 
Robinson bound (11). 

We envisage the bound (11) to be particularly useful for hnite systems of 
intermediate size where the matrix exponential can be computed numerically with 
reasonable effort. However, since (11) is sharper than the bounds in [23, 25], a 
thermodynamic limit will exist (at least) under the same conditions required in those 
proofs, and in particular for H-dimensional regular lattices with power law interactions 
and exponents a > D. 


Appendix A.2. Lieb-Robinson bounds in rescaled time 

As in section 3, we consider a H-dimensional lattice A consisting of N sites, a Hilbert 
space (2) with hnite-dimensional local Hilbert spaces, and a generic Hamiltonian H with 
n-body interactions (3). We require that H satishes conditions (4) and (5). For the 
proof of a Lieb-Robinson-type bound, we follow the general strategy of [25], augmented 
with the tHA-rescaling taken from [30]. 

As a shorthand we introduce 

X{t) = [OA{t),OB]. (A.5) 

for the commutator on the left-hand side of (8). Differentiating with respect to t yields 
x'it) = i [x{t), lA{t)] + i pAit), [lA{t), Ob]] , (A.6) 
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where I a = is the set of local Hamiltonian terms that have non-zero 

overlap with A. Using the boundedness of OA{t) we apply Lemma A.l of Ref. [25] to 
the norm-preserving first term of (A.6), yielding 


llx(i)l|-||x(0)||<2||O4|| 5^ 

z-.zr\A^% • 


A\ 


\\[hzis),OB]\\ds. 


Next we define 


Cos{A-t)-= sup 

Oa&Oa IlfdAll 


(A.7) 


(A.8) 


where Oa is the set of observables compactly supported on A. Making use of this 
definition, (A.7) can be rewritten as 


Cob ~ ddog ( 0 ) t) 


^ E 




\\hz\\CoB{Z,s)ds. 


Eq. (A.9) can be applied recursively to show that 

Co.(Ai)< 2 IIOBllf: 

with coefficients 


n=l 


nl 


— 


X] ••• n \\^zi\\5B{Zn), 


(A.9) 


(A.IO) 


(A.ll) 


ZlCA Z2CA Z„CA l=l 

ZinA7t0 ZjnZjT^a z„nz„_i7^0 


Sb{Z) = 


where 

'0 if Z n R ^ 0, 

1 otherwise. 

Under the conditions (4) and (5) these coefficients can be bounded by 

pn-i^n 

Inserting (A. 13) into (A.IO) and using the definition (9) of rescaled time r, one obtains 

2 ||Ob|||A||R| 


(A.12) 


(A.13) 


Cos(A,f) < 
and this implies the bound 


p{i + d{A,B)y 


(exp [2pA |r|] - 1), 


\\[Oa{t^k),Ob\\\< 


2 || 0 ^|||| 0 s|||A||R| 


(exp [2pA |r|] - 1) 


p{l + d{A,B)y 

in rescaled time r, valid for power law interactions with exponents a > 0. 


(A. 14) 


(A.15) 
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Appendix A.3. Discussion of the bound in Reference [26] 

In [26] a Lieb-Robinson-type bound was derived whose functional form consists of a 
linear (cone-like) and a faster-than-linear (supersonic) contribution. This bound is a 
major improvement over that in [23] in the regime of large a, where the former becomes 
more and more similar to a nearest-neighbour bound, as it should. Here we scrutinize 
the applicability of the bound in [26] for describing the cone-like and supersonic features 
of long-range models with intermediate exponents a (roughly in the range 3 < a < 8). 
The bound in [26] is derived for Hamiltonians 


with two-body interactions hij satisfying 

Hop'll < 


(A.16) 


dii^y 


(A.17) 


on H-dimensional regular cubic lattices. For exponents a > 1 a bound of the form 

\mt),B]\\ 


< Ti + Ta 


(A.18) 


2|A| \B\ 

is obtained, where A and B are observables on lattice sites that are a distance 6 apart, 
and 


Ti = ci- 


^Vit 


3/^(5 


T2 = C21 


^V 2 t 


with Cl = A“^, vi = 2A^e, C 2 = (A9^)“^, V 2 = 2A^9'^, A = and 0 < p < 1. 

Ti has the same functional form as the classic Lieb-Robinson bound for Hamiltonians 
with hnite-range interactions [22], which is known to produce a linear, cone-shaped 
causal region. T 2 has the functional form of the bound originally derived by Hastings 
and Koma [23]. Both, Ti and T 2 contain the free parameter p, which determines, among 
other things, the slope of the linear soundcone. So the “velocity” associated with the 
cone can be tuned to an arbitrary value, irrespectively of the physical behavior of the 
model studied. 

Based on (A.18) and (A. 19), the sharpest bound 

/ glJlt _ ^ gD2t _ ^ \ 

B(i.i) = ™u(c,^^ + c,[(i_^)^l„) (A.20) 

is obtained by minimizing, separately for each value of 6 and t, over the free parameter 
p. From the contour plots of B in hgure Al it becomes clear that the “linearity” of Ti 
can be deceiving, as a linear, cone-like regime is not particularly prominent, not even for 
larger exponents like a = 8. Of course it is always possible to construct a linear-looking 
bound by weakening R, but this would be unr elated to the physical behavior of the class 
of models studied. 

Another, more sophisticated bound has recently been put forward in [27], but the 
form of the propagation front has not yet been analyzed and discussed (beyond the 
long-distance asymptotics). 
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Figure Al. Contour plots of the bound (A.20) in one spatial dimension for a = 6/5 
(left), 4 (center), and 8 (right). Even for a = 8 there is at best a hint of a linear 
regime. See the matrix exponential bound in figure 3 for comparison. 


Appendix B. Information capacity of the long-range Ising model 

In this appendix we prove that supersonic transmission through a quantum channel can 
occur for any a < 2, also for measurements performed on single lattice sites. Like for 
the study of the group velocity of the long-range hopping model in section 4.3, we find a 
threshold value of a = 2 below which propagation becomes supersonic. The proof uses 
techniques from Ref. [15] and applies them to a slightly more involved model for which 
supersonic propagation is found to occur also for single-site measurements. 

We consider a finite one-dimensional lattice A = {1,..., A} consisting of N sites. 
To implement a quantum channel, we encode a signal on site 1, and measure the effect of 
that encoding after a time time at site N. On this lattice we define an Ising Hamiltonian 
with arbitrary couplings, 

^ = E (■3.1) 

i<j 

Defining the sublattices A = {1}, B = {N}, and S = A\{AU B), the Hamiltonian can 
be rewritten as 

H = Has + Hab + Hsb + Hss (B-2) 

with 

-EE (B.3) 

iex jeY 

where X,Y E {A, S,B}. As an initial state we choose 

p(0) = I Dl idl I 'Ds sdl ® |+)Af Af(+| (B-4) 

seS 

with af |t)i = It)*, < lt)i = - It)*, and |+)^. = (|t)j + |t)j)/V2- Initially all the spins 
are pointing down, except the one ai B = {A}. 

A binary quantum channel is implemented by starting the time evolution either 
with p(0) (sending a “0”), or starting with UApiO)U\ (sending a “1”), where Ua is a 
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unitary supported on A only. The classical information capacity Ct can be bounded 
from below by the probability to detect, by measuring according to a positive operator 
valued measure ttb, a signal at B after a time t, 

Ct>Pt = |Tr {Nt [p(0)] ttb} - Tr {Tt [p(0)] tib]1 (B.5) 


with 


Nt [p(0)] := TrA\B e ' p(0)e 




Tt [p(0)] := TrA\s 




(B.6) 

(B.7) 


In the following we compute a lower bound on the right-hand side of (B.5), and study 
this bound as a function of the channel length, i.e., the distance between A and B. 

We choose ttb = |+)Ar iv(+| and Ua = | li where the latter is a spin flip 

operator on the first lattice site. For the time-evolved density operator in (B.6) we hnd 

(B.8) 




pit) = e 

p V-1 
X 

- S=1 

N-l 


4-)ssU I ® \+)n w(+| 




I exp 


S=1 


■ N-l 

JrN(^N 

r=l 


|+)Ar Ar(+| exp 


Af -1 


JrN^^ ]\f 


r=l 


All the exponentials not supported on B add up to zero since the initial state prepared 
on A\B is an eigenstate of the Ising Hamiltonian. Taking the trace gives 


Tr {Nt [p(0)] TTs} = - <{ 1 -h cos 


2 t E JrN + JlN 


vrSS 


A similar calculation shows that 


Tr [Tt [p(0)] TTs} = - <1 1 + cos 


2t ( ^ JrN 


J 


IN 


ktGS 


The probability of detecting a signal in B at some time f > 0 is then given by 

cos 2 t I 
L \res 


Pt = 


2 t E JrN + JlN 


— COS 


2 t E JrN — JlN 


■^r£S 


sm 


2 t ^ JrN I 

^ r-es / 


sin (2f Jiat) 


(B.9) 

(B.IO) 

(B.ll) 

(B.12) 


To derive a nontrivial (nonzero) lower bound on pt, we target the regime before 
oscillatory behavior in (B.12) sets in. Using the inequality 


‘2x 

sin X > — for 0 < a; < 7r/2 

TT 

and assuming power law interactions Jtj = |f — j|~“, we obtain 


Pt > 


At 


4t 


N-l 

E 


71 [N — 1)" 71 [N — r)“ ’ 


(B.13) 


(B.14) 
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valid for times 


N-2 . 
TT / V—\ 1 

^ < t/ / — 
4 / r“ 

r=l 


Interpreting the snm in (B.14) as an npper Riemann sum, we have 


N-l 

E 


Then we can bound pt by 
Pt > 


N-2 

E 

r=l 


-V-2 


- > 


dr 


(r + !)“■ 


7r2(a — 1) (A^ — 1)^ 


=-Pf 


(B,16) 


(B,16) 


(B.17) 


For a > 1 and large N the second term in the square bracket in (B.17) is much smaller 
than 1, and we obtain 

let^ 1 

-t ^ 7r2(a-l) (iV-l)“ ^ 

for the large-iV asymptotic behavior of the bound p^. In our setting, d = N — 1 is the 
distance between the regions A and B. To determine the shape of a contour line at 
which p^ is equal to some constant e, we set 


^ = Pt^Jc 


and we can read off that 
5 oc 


(B.19) 


(B.20) 


along any of those contour lines. Eq. (B.20) describes faster-than-linear (supersonic) 
growth of d for a < 2. It is straightforward to extend the above calculation to more 
general initial conditions as well as to lattices of arbitrary dimension. 
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